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Abstract. Two prominent magnetoelectrics MnWCU and CuO possess low- 
temperature commensurate paraelectric magnetically ordered phase. Here using Monte 
Carlo simulations we show that the walls between the domains of this phase are 
ferroelectric with the same electric polarization direction and value as those in the 
magnetoelectric phases of these compounds. We also suggest that experimental 
observation of electric polarization of domain walls in MnWCU should help to determine 
the macroscopic interactions responsible for its magnetoelectric properties. 
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1. Introduction 

Increasing demand in miniaturization of functional elements of devices has tremendously 
stimulated research in the physics of domain walls in recent decades |T|. The walls 
separating different domains are often found to have new properties different from 
those of the domains. The examples of such exploitable phenomena are ferroelectric 
domain walls separating paraelectric domains 13 El (a u i, conducting domain walls 
between insulating domains [7], and superconducting twin walls [8]. Thus, domain walls 
provide new physics, which can be used to create new types of devices, such as, for 
example, the racetrack memory [9j and magnetic domain wall logic [10.] or can be used 
as injectors/detectors of spin currents [IT] . 

At the same time the last decade has seen rapid increase of interest in materials 
showing multiferroic properties ra- Ferroelectromagnets (i.e. materials showing 
independent magnetic and ferroelectric orders) and magnetoelectrics (i.e. materials 
displaying magnetically induced ferroelectric order) are of particular interest due to 
the possibility of creating new types of devices such as, for example, devices for 
reading, writing, and storage of information [[13], various sensors [14], microwave [15] 
and spintronics [16] devices, and wireless energy transfer and energy harvesting 
technologies [15] . 

Technological applications of multiferroics discovered to date, however, are still 
significantly hindered by a number of their properties. Ferroelectromagnets usually 
show high temperatures of ferroelectric and magnetic phase transitions and high electric 
polarization. However, the usually large difference between the temperatures of the 
magnetic and ferroelectric phase transitions and the difference in the sources of the 
two orders result in low magnetoelectric response ra- In turn, magnetoelectrics show 
direct coupling between magnetic order and induced electric polarization, but possess 
low magnetic phase transition temperatures and low electric polarization values. The 
magnetoelectric with the highest phase transition temperature discovered to date is 
cupric oxide CuO [17], which displays modulated magnetic order and ferroelectricity 
between approximately 210 and 230 K. 

The microscopic nature of the magnetoelectric effect is usually interpreted using the 
inverse Dzyaloshinskii-Moriya interaction [18] or the spin-current model [19]. In both 
models the electric polarization arises as the result of interaction of two neighboring 
noncollinear spins. However, single ion contribution to the microscopic origin of the 
magnetically induced polarization was shown to exist [20]. 

Common occurrence of modulated magnetic structures in magnetoelectrics has led 
to pervasive notion that incommensurate magnetic order such as, for example, spiral 
or cycloidal is necessary to induce electric polarization [2T]J. However, ferroelectricity 
can be induced by commensurate magnetic structures, which is observed, for example, 
in RM112O5 (R = rare-earth) magnetoelectrics [22] . The intimate coupling between 
modulated magnetic order and magnetically induced electric polarization can be 
understood as the symmetry property of the relevant magnetic order parameters, which 
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allow Lifshitz invariants responsible for modulation and induce improper ferroelectric 
phases [ 23] . 

Two approaches exist to macroscopically describe the magnetoelectric interaction. 
The first one consists in inclusion of the terms in the thermodynamic potential, which 
can be schematically written as 
drj 2 


PaVl 


dp' 


( 1 ) 


where a and ft are coordinates, P a is the electric polarization component, and r)\ and 
r /2 are either different components of a single magnetic order parameter or two different 
magnetic order parameters. Such magnetoelectric effect due to inhomogeneous magnetic 
order was first studied in [23], where the calculation of electric polarization of Bloch and 
Neel magnetic domain walls was outlined. Subsequently, the same ideas were used to 
describe the bulk magnetoelectric effect in spiral magnets [25] and led to the appearance 
of the term flexomagnetoelectric effect [26], 

The second approach consists in considering the contribution to the thermodynamic 
potential of the form 


P a r]ir]2. 


( 2 ) 


Such macroscopic coupling was successfully used to interpret magnetoelectric 
phenomena in a number of magnetoelectrics [27], [28, '291 ED]. Apparently 
both interactions © and (J 2 [) can describe magnetoelectric phenomena in many 
magnetoelectrics. However there is an important difference between them consisting in 
the fact that ([I]) requires spatially varying magnetic order parameters to induce electric 
polarization. The experimental fact that ferroelectricity occurs mostly in spatially 
modulated magnetic structures together with the growing number of crystals showing 
polar magnetic domain walls, which is explained using flexomagnetoelectric interaction, 
give evidence in support of the interaction (0Q) . However, the interaction ([2]) is capable of 
describing magnetoelectric effect in both incommensurate and commensurate magnetic 
structures and, therefore, provides equal footing for such description. 

The purpose of this work is twofold. Firstly, we show that the walls separating the 
domains of low-temperature commensurate paraelectric magnetic phases of MnWCU 
and CuO are polar. Secondly, we show that this fact can be used to experimentally 
distinguish between the cases of interactions ([T]) and ([ 2 ]) in theoretical description of 
magnetoelectric properties of M 11 WO 4 . The paper is organized as follows. In section [2] 
we describe the spin Hamiltonian and the results of the Monte Carlo simulations, in 
section [3] we discuss the obtained results, whereas in section [4] we present the conclusions 
of our work. 


2. Spin Hamiltonian and Monte Carlo studies 

M 11 WO 4 possesses a monoclinic structure with the monoclinic angle [3 ~ 91° described 
by the space group P2/c (C^J [31]. Upon lowering the temperature it undergoes 
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Table 1 . Magnetic superexchange constants for the Mn 2+ ions in MnW 04 given in 
units of meV. The Mn - Mn distance is given in A. 

J\ J 2 J 3 J& J 5 


Mn - Mn 3.28 4.4 4.82 4.99 5.75 5.8 

|34| -0.168 -0.116 -0.364 0.356 0.018 -0.438 

[33 -0.84 -0.08 -0.64 -0.52 0.1 -0.86 



J7 

A 

J 9 

J 10 

^11 

D 

Mn - Mn 

5.87 

6.49 

6.56 

6.88 

7.01 


m 

0.02 

0.424 

-1.96 



-0.122 

ESI 

-0.24 

0.04 

-0.52 

-0.3 

0.04 

-0.18 


a sequence of magnetic phase transitions at T N =13.5 K, T 2 =12.7 K, and 7\=7.6 K, 
which lead to the appearance of magnetically ordered states AF3, AF2, and AF1, 
respectively [32]. The structure of the low-temperature commensurate magnetic phase 
AF1 is described by the wave vector k c = (|, |), whereas the incommensurate phases 

AF2 and AF3 by the wave vector k inc = (—0.214,|,0.457) [3T] . The phase AF2 is 
ferroelectric with polarization along the b axis [32]. In the phases AF1 and AF3 without 
applied external magnetic fields the spins in M11WO4 are directed along the easy axis 
in the ac plane making an angle of about 34° with the a axis [33]. I 11 the phase AF2 
additional spin component along b appears. In the following we assume that the axes 
x, y, and z are directed along the easy axis, b axis, and perpendicular to the easy and 
b axes, respectively. 

For the analysis of the magnetic structure of domain walls in M 11 WO 4 we employ 
Monte Carlo simulations of the Heisenberg Hamiltonian 

H = ~\ J{R + t,R'+ t')S(R + t)S(R'+ t') 

Z R,R',t,t' 

+-^2 D a Sa(R + t), (3) 

R,t,a 

where S(R + t) are classical spins of length |, R and R' are lattice vectors, t, t' = ti, t 2 
give the positions of the two Mn 2+ ions in the unit cell, J(R + t, R' + t') is the 
superexchange constant between the spins at R + t and R' + t', D a is the single-ion 
anisotropy, and a — x,y, z. The first sum in (J3]) is evaluated over pairs of spins. 

Inelastic neutron scattering was used to determine up to eleven superexchange 
constants in M 11 WO 4 [34, [35] . which are summarized in table [I] The free energy of 
M 11 WO 4 corresponding to the Hamiltonian ([3]) and obtained using mean field approach 
was thoroughly studied using Fourier transform and was found to give a correct 
phase transition sequence and good correspondence with the experimental data [36]. 
Furthermore, both constants sets of table |T] are shown to give results in qualitative 
agreement with the experimental data on M 11 WO 4 [36]. For our Monte Carlo simulations 
below we use the exchange constants of [35] with D x = D, D y = 0, and D z = —D. 
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The Monte Carlo studies were performed using the Metropolis algorithm. At each 
temperature the system was allowed to relax for 5-10 9 Monte Carlo steps, after which 
the data was averaged for 50-10 6 steps. This approximately amounted to 8.7-10 5 and 
8.7T0 3 steps per spin, respectively. 

The phase transition sequence and the appearing magnetically ordered states in 
M 11 WO 4 can be described by magnetic order parameters belonging to the wave vector 
k c [37]. Furthermore it was argued that the praphase approach to the description of 
phase transitions in MnW 04 has significant advantages [30]. Thus in the following we 
use magnetic order parameters belonging to the wave vector k c of the orthorhombic 
praphase with the space group Pmcm (D^) [30]. The hypothetical phase transition 
Pmcm - P 2 /c is described by the component Uxz of the homogenous deformation 
tensor, which has to be assumed nonzero value for the description of the monoclinic 
phase. In the following the axes X, Y, and Z are assumed to be directed along the a 
axis, the b axis, and perpendicular to the a and the b axes, respectively. Choosing the 
praphase as the starting point for determination of the magnetic order parameters has 
no particular influence on the obtained results, but makes the present study consistent 
with our previous work on M 11 WO 4 . 

In the orthorhombic praphase structure M 11 WO 4 possesses two Mn 2+ ions in the 
unit cell located at (|, 0.6853, -j-) and (|, 0.3147, |), whose spins are denoted by Si 
and S' 2 , respectively. The magnetic representation in k c point of the Brillouin zone 
is given by 3Pi, where Pi is a four-dimensional irreducible representation |30j . Upon 
the hypothetical phase transition Pmcm - P 2 /c Uxz splits Pi into two irreducible 
representations, which are denoted by Gi and G 2 in the monoclinic structure. Denoting 
by (gia, 92a, 93a, g±a) (a = x, y, z ) the magnetic order parameters transforming according 
to Pi and describing the spin components along the respective directions we obtain the 
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manganese spins in a unit cell given by nidi + n 2 a 2 + n 3 a 3 as 


Sly = -Z 


Six — — £ (dlx + 92x — 9Ax — 9Ax) COS j 

. fUiTl\ 

— \9lx — 9'2 x — 9Ax + 9 Ax) Sill ^ —— J 
i.9ly ~ 92y ~ 93y + 9Ay) COS ^ ^ ) 

+ {.9ly + 92y ~ 9Ay ~ 9Ay) sin 

Su = £ ( 9iz + 92 z + 9az + 9az ) cos — j 

-(giz - 92z + g 3z - 9Az) sin (^y) 
S'2x £ ( 9lx 92x T 9 3 x 9Ax) COS ^ ^ 

+ (9ix + , 92 .x + 9Ax + gAx) sin ) 

S r 2y = — ^ (9ly + 92y + 9Ay + 94?;) COS — j 
~(9ly — 92y + 9Ay ~ 9Ay) sin ^ ^ ) 
Szz = f (gu - 92 z - gAz + gAz) COS 

+ (9ix + 92 z - 9az - 9az) sin 


(4) 


where ai, a 2 , and a 3 are Bravais translations and f = (—i)™ 2 + n3 . 

Previous analysis of the exchange symmetry in M 11 WO 4 showed that the 
set of order parameters (gi a , 92 a, 93 a, 9 ao ) splits into two exchange multiplets 
{93x, gAx, giy, g 2 y, g 3 z, 9Az) and {gix,g2x,g3y,9Ay,giz,g2z) [30]. In M 11 WO 4 the magnetic 
structures without applied magnetic field belong to the first exchange multiplet [31] 30]. 
Further splitting within exchange multiplets consistent with the crystallographic 
symmetry is due to single-ion anisotropy D a . 

The theoretical model of phase transitions in M 11 WO 4 [30j, which is in accordance 
with experimental data [33113T1138] . reveals that the magnetically ordered phases AF3 
and AF2 are described by spatially modulated order parameters (g 3x (r), g^xif)) an d 
(9iy(r), 92 y(r), 93 x(r), gAx(r)), respectively, whereas the commensurate phase AF1 by 
the order parameter (g 3 x ,gA X ) = ( 9 , 0 ). Therefore, the phase AF1 can be described by 
four domains with ( g 3x , g ix ) equal to (±9, 0) and (0, ±9), which are shown in Figs. Oja) 
and (b). In this work we study the domain walls that separate these domains and are 
oriented perpendicular to the crystal axes a and c, which we hereafter refer to as a-walls 
and c-walls, respectively. From the symmetry of the problem there are essentially three 
kinds of walls, those that separate the domains (9,0) and (0,9), the domains (9,0) and 
(—9, 0), and the domains (9, 0) and (0, —9), which we hereafter denote type I, II, and III 
domain walls, respectively. 

Firstly, however, we study the sequence of magnetic phase transitions in M11WO4 
and for simplicity consider only the spatial variation of magnetic order along the a 
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Figure 1 . Ordering patterns of Mn 2+ spins for various values of the order parameter 
(gi y , g 2 V , gzxi 9 ax)- The ordering pattern for negative value of g is obtained from that 
for the positive value by spin inversion. 


axis by performing Monte Carlo simulations of a system with dimensions 80 x 6 x 6 
unit cells and periodic boundary conditions. Figure EJa) shows spatial variation of 
the order parameter (gi y , g 2 y , g^x-, 9 a x ) in a decreasing temperature run. ft can be 
seen that upon decreasing temperature a magnetic structure with spatially modulated 
order parameter (0, 0, g 3x (X), g ix (X)) appears at T N . We associate this phase with 
AF3. The modulation period is approximately 20 unit cells giving the wave vector 
component of 0.2 along a, which roughly coincides with the experimental value 0.214. 
Simulations of larger systems (1000 x 6 x 6) gave approximately the same modulation 
period. The phase AF3 is followed by the appearance of the y component of spins 
(giy(X), g 2y (X), g 3x (X), g 4x (X)) and the resulting structure is associated with the 
phase AF2. Further temperature decrease does not result in a phase transition to 
a commensurate phase. We attribute the absence of such transition in our Monte 
Carlo simulations to its strongly first order character and suggest that further studies 
of this phase transition are needed, which is not in the scope of the present paper. 
However the commensurate state AF1 is stable when chosen as an initial state. The 
respective simulations of the commensurate state during increasing temperature run 
are shown in figure [2](c). It can be seen that the value of the nonzero component 
of the order parameter g 3x gradually decreases with increasing temperature and at 
temperature between approximately 21.7 K and 24.6 K a first order phase transition to 
an incommensurate phase occurs. According to our simulations at T = 0 K the internal 
energy of the commensurate state is lower than that of the modulated state, which 
ensures a phase transition to the commensurate state at temperature T\ intermediate 
between T = 0 K and T 2 . 

The macroscopic symmetry allows magnetoelectric interactions responsible for the 
emergence of electric polarization along the crystal b axis, the lowest of which in powers 
of magnetic order parameters and spatial derivatives can be written in the form 


*b\9lx93y + 92x9 Ay + 9 3: 
( dg 4y dg 3y 
w' 7 OX 92x 8X 




, „ 9<hy \ 
+ 93x gx 9 ax gx J , 
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Figure 2. (a) Spatial dependence of (gi y , g 2y , g 3 x, 9 Ax) along the a axis at decreasing 
temperature run. (b) Spatial dependence of A, I 2 , and I 3 at different temperatures 
corresponding to Figure (a), (c) Spatial dependence of (ffiy, g 2y , g 3 x, gr x ) along the a 
axis during increasing temperature run, when AF1 is chosen as the initial state. The 
horizontal axes in all figures are in cell units. 
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dX dX dX 
Invariants similar to (} 6 j) — ([9}) exist, in which every term in the parantheses of the form 
gixdgjy/da should be substituted by gj y dg ix /da, where i , j = 1,2, 3,4 and a = X,Z. 
According to our simulations and experimental data the relevant order parameter in 
M 11 WO 4 is ( 9 i y , g- 2 y, g 3x , 941 ), with all other components of the order parameters { 9 * 7 } 
(i = 1,2,3,4, 7 = x,y,z) experiencing only statistical fluctuations. Therefore, the 
magnetoelectric invariants (l5]l — ()9|) can be written in simpler forms, which give electric 
polarization P b proportional to the following combinations of order parameters 


p b ~ 

(93x9ly T 94x’92y) 

— h, 

Pb ~ 

( dg 2 y 

V /3x da 

n d9l A 

94x da ) 

= h, 

Pb ~ 

( dg 3x 

^9 Ax \ 

= 7.3, 

K da 

9ly da ) 


( 10 ) 

( 11 ) 

( 12 ) 


where a = X, Z. Figure EJ^b) gives I\, I 2 , and I 3 for several temperatures as functions of 
the cell number along the a axis. The dependencies experience a sawtooth character and, 
therefore, the electric polarization should be proportional to average values, which are 
indicated by solid lines. Depending on the phase shift between (gi y ,g 2y ) and ( g 3x ,gA X ), 
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Figure 3. (a) Temperature dependence of amplitudes of g\ v , g 2 V , gz x , and g^x 

in incommensurate states during decreasing temperature run. (b) Temperature 
dependence of the average values of I±, I 2 , and I 3 during decreasing temperature 
simulation. 

which is achieved during cooling of the system, the quantities / 1 , J 2 , and I 3 can have 
both positive and negative values giving both directions of electric polarization 

Figure E^a) shows temperature dependencies of the amplitudes of modulated order 
parameter components g± y , g 2y , <? 3x , and g± x during decreasing temperature run. From 
these data one can estimate the phase transition temperatures as TV ~ 26.7 K and 
T 2 fa 23.2 K. Our Monte Carlo studies give transition temperatures higher than those 
observed in experiments. We attribute this discrepancy to inaccuracy in exchange 
constants, which, however, from our point of view does not affect qualitative results 
obtained in the present work. Figure |3|b) gives the temperature dependence of the 
average densities of Ji, J 2 , and J 3 , which are essentially proportional to the electric 
polarization. In agreement with experimental data P& appears at T 2 when all components 
of the order parameter (gi y , g 2y , g 3 X , 9 ix) condense. 

Next, we study the walls between the domains of the commensurate paraelectric 
phase AF1. Figure 0] shows variations of the order parameter and polarization across 
the domain walls of types I, II, and III, which are oriented perpendicular to the a axis 
of the crystal. For this purpose a system with dimensions 80 x 6 x 6 unit cells was 
used. It can be seen that all domain wall types are ferroelectric. The values of 1 1 , J 2 , 
and I 3 inside the walls are of the order of those in the incommensurate phase AF2 
shown in figure E|b), which allows concluding that the value of polarization of the walls 
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Figure 4. (a) Spatial variation of (gi y , giy, g 3 x, 9ix) across the type I a-wall. (b) 
Spatial variation of I±, Ii, and I 3 across the type I a-wall. (c) and (d) are same as 
(a) and (b), but for the type II wall, whereas (e) and (f) are for the type III wall, 
respectively. The horizontal axes in all figures are in cell units. 


should coincide with that of the bulk value in the phase AF2. Similar to the case above, 
depending on the signs of g ly and g-i y the quantities E, I 2 , and I 3 can have both positive 
and negative values, which can result in both directions of electric polarization P,\ of the 
wall. The states with positive and negative Pb have the same energy. 

The order parameter and electric polarization profiles across domain walls oriented 
perpendicular to the crystal c axis are qualitatively similar to those of figure [4] and will 
not be given here. Figure [5fa) shows the temperature dependence of the domain wall 
thicknesses for both wall orientations, which we calculate as full width at half maximum 
of electric polarization. It can be found that the thickness of the a-walls, which is of 
the order of 7 - 15 unit cells, is larger than than that of the c-walls, for which it takes 
values of the order of 5 unit cells. With increasing temperature the thickness of the 
walls gradually increases, whereas average polarization, which is shown in figures Efb) 
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Figure 5. (a) Temperature dependence of domain wall thickness. Empty and filled 
signs show data for domain walls oriented perpendicular to the crystal axes a and 
c, whereas triangles, squares, and circles represent domain wall types I, II, and III, 
respectively, (b) and (c) give temperature dependence of average values of I \, E , and 
I 3 for domain walls oriented perpendicular to the crystal axes a and c, respectively, 
whereas filled triangles, empty squares, and empty circles denote data for domain wall 
types I, II, and III, respectively. 


and (c), gradually decreases. The temperature evolution of the properties of domain 
walls in figure [5] are shown for temperatures up to approximately 17 K, because at 
higher temperatures the walls become less discernible and a phase transition to an 
incommensurate state occurs. 

3. Discussion 

In the previous section we have shown that the walls between the domains of the 
paraelectric phase AF1 of MnW 04 are ferroelectric with electric polarization along the 
b axis of the crystal, which coincides with its direction in the phase AF2. It is known 
from the previous models of phase transitions in M 11 WO 4 that I\ can only appear when 
both parts (, gi y ,g 2y ) and (g 3x ,g 4 x) of the order parameter (g ly , g 2y ,g 3x , g*x) condense, 
which follows from the form of magnetoelectric interactions (l5]l — ()9|) | 30l [28 ]. In our 
simulations we find that inside the walls both parts of the order parameter condense 
in such a way that the wall becomes polar. The appearance of nonzero (gi y , g 2 y) inside 
the wall can be explained in the following simple way. Consider, for example the type I 
wall. The change from one domain with the order parameter (0, 0,g,0) to the other 
with (0,0, 0,(7) without nonzero ( gi y ,g 2 y ) would require passing through a state with 
the order parameter (0,0 ,g',g'). Such state, however, results in ordering of only one 
of the two spins Si and S 2 in every unit crystal cell as follows from (J4J) and is shown 
in Fig. (He). (Figure QJd) shows the phase state (0,0 ,g',—g')-) Therefore, this state 
possesses considerably higher exchange energy and is less favorable than the state with 
nonzero (gi y ,g 2y ) provided that the single-ion anisotropy is sufficiently small, which is 
the case in M 11 WO 4 . In terms of irreducible representations of the monoclinic space 
group the phase AF1 is described by G 2 , whereas Gi additionally appears inside the 
walls between the domains of the phase AF1. Simultaneous appearance of Gi and 
G 2 inside the wall results in its ferroelectric polarization ffl. The ordering patterns of 
Mn 2+ spins in the middle of the type I, II, and III walls, which approximately correspond 
to the phase states (- g , - g , - g , -g), ( 0 , -g, 0 , - g ), and (g, -g , g , - g ), respectively, are 
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Figure 6. Ordering patterns of Mn 2+ spins in the middle of the type I (a), type II (b), 
and type III (c) walls. The figure was created using the VESTA software [59] , 


shown in Fig. [ 6 j 

We have also identified three possible macroscopic magnetoelectric interactions, 
which give electric polarization proportional to Ji, J 2 , and J 3 . It follows from 
our Monte Carlo studies that all three possible magnetoelectric interactions can 
contribute to ferroelectric polarization in the phase AF2 as well as in domain walls. 
However, it could be the case that one or two of them prevail over the others. 
Besides, these interactions give different results for commensurate states, i.e. the 
interactions possessing spatial derivatives (J 2 and J 3 ) cannot give electric polarization 
in commensurate magnetically ordered state, whereas I\ can induce polarization in 
commensurate as well as incommensurate magnetic states. Therefore, it is essential to 
determine which of these interactions are responsible for magnetoelectric properties of 
M 11 WO 4 . Here we argue that experimental studies of the walls between the domains of 
the phase AF1 can clarify this problem. Indeed, the crystal structure of M 11 WO 4 
possesses ...-Mn-O-Mn-O-... zigzag chains running along the c axis, whereas along 
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the a and b axes the interaction between manganese spins is more complex (Mn-O- 
W-O-Mn). Thus, the exchange paths along the a and the c axes are substantially 
different. Therefore, one can argue that if the magnetically induced electric polarization 
in M 11 WO 4 is due to differences of the order parameter in neighboring cells, i.e. if 
the relevant magnetoelectric interactions in M 11 WO 4 are proportional to the spatial 
derivatives of the order parameter as in J 2 and J 3 , then the a-walls should possess 
zero or low electric polarization in contrast to the c-walls, which should be ferroelectric 
with electric polarization comparable to that in the bulk. On the contrary, if similar 
polarization is revealed experimentally for both the a- and the c-walls, then one can 
argue that the relevant magnetoelectric interaction is ([5]). 

In the available literature the studies of the domains or domain walls in M 11 WO 4 
are limited to those in the ferroelectric phase AF2 and are mostly concerned with the 
dynamics of electric polarization switching by electric or magnetic fields dmmuz]. 
In [43] the effect of poling electric field on ferroelectric properties of M 11 WO 4 was 
studied. The authors observed nonzero electric polarization below T\ when the sample 
was cooled down from T T N under poling electric field of 10 kV cm -1 to temperatures 
below Ti, after which the field was removed and a residual polarization was observed. 
This phenomenon was interpreted as enhancement of stability of the phase AF2 by 
electric field, which results in surviving of the clusters of ferroelectric phase AF2 in 
the paraelectric phase AF1. The authors estimated the volume fraction of the phase 
AF2 at T = 4 K to be as high as 13% by comparing the polarization value to that 
of the pure phase AF2. We agree with the authors of |43| that electric field enhances 
the stability of the phase AF2. This enhancement, however, is present only when the 
poling field is applied to the sample and should disappear upon its removal, which at 
T < Ti should result in a phase transition of the remains of the phase AF2 to the 
phase AF1 and disappearance of polarization. From our point of view the presence of 
electric polarization at temperatures below Ti can be interpreted as polarization of walls 
between the domains of the phase AF1. The poling electric field present upon the phase 
transition to paraelectric phase AF1 from AF2 should (i) pole the appearing domain 
walls and align their polarization along the field, and (ii) should in general increase the 
number of domains and, respectively, the number of the walls between them. Therefore, 
we argue that these experimental results can be interpreted as supporting the picture 
of polar domain walls in the phase AF1 suggested in the present work. 

The analogy in the phase transition sequences of M 11 WO 4 and CuO, which is also 
reflected in similarity of their macroscopic descriptions [30] 44], argues that the same 
phenomena suggested for the low-temperature phase of M 11 WO 4 in the present work 
can take place in the phase AF1 of CuO. Electric polarization of domain walls in the 
phase AF1 of CuO will extend the exploitable temperature range of magnetoelectric 
properties of CuO down to 0 K from the rather narrow range 213 - 229 K of stability 
of the phase AF2. 
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We have shown that the walls between the domains of the low-temperature 
commensurate paraelectric phase AF1 of M 11 WO 4 and CuO are ferroelectric. The 
value and direction of electric polarization inside the walls coincide with those of the 
magnetoelectric phase AF2 of the respective compounds. Experimental observation of 
electric polarization of domain walls in M 11 WO 4 should help clarifying the macroscopic 
interaction responsible for its magnetoelectric properties: flexomagnetoelectric or 
trilinear. Electric polarization of domain walls of CuO significantly extends the 
temperature range of its exploitable magnetoelectric properties. Together with the 
possibility of switching the electric polarization of the walls and the complex behavior 
of these magnetoelectrics in magnetic fields this may open new interesting phenomena 
in the physics of these magnetoelectrics. 
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